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Abstract. On the basis of a semi-classical analysis of vacuum energy in an 
expanding spacetime, we describe a non-singular cosmological model in which 
the vacuum density decays with time, with a concomitant production of matter. 
During an infinitely long period we have an empty, inflationary universe, with 
H Ri 1. This primordial era ends in a fast phase transition, during which H and 
A decrease to nearly zero in a few Planck times, with release of a huge amount 
of radiation. The late-time scenario is similar to the standard model, with the 
radiation phase followed by a long dust era, which tends asymptotically to a 
de Sitter universe, with vacuum dominating again. An analysis of the redshift- 
distance relation for supernovas la leads to cosmological parameters in agreement 
with other current estimations. 



1. Introduction 

The role of vacuum energy in cosmology has acquired a renewed importance, with 
current observations pointing to the existence of a negative-pressure component in 
the cosmic fluid. This has reinforced the cosmological constant problem pQ, whose 
origin can be understood if we write the energy density associated to vacuum quantum 
fluctuations. In the case of a massless scalar field, it is given by the divergent integral 



Jo 

This integral can be regularized by imposing a cutoff m in the superior limit of 
integration, leading to Ao ~ to 4 . The same result can be derived by introducing a 
bosonic distribution function in (fl]), 



which is equivalent to consider the vacuum fluctuations thermally distributed, at a 
characteristic energy m. 

For a cutoff of the order of the Planck mass, this leads to a vacuum density 
120 orders of magnitude above the currently observed dark energy density. One may 
argue that m should in fact be much smaller, because vacuum fluctuations above the 
energy scale of the QCD phase transition - the latest cosmological vacuum transition 
- would lead to quark deconfincmcnt. However, even with this value for m, we obtain 
a vacuum density 40 orders of magnitude above the observed value. 
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The cosmological constant problem may be alleviated by the following reasoning. 
The above energy density is obtained in a flat spacetime, but in such a background 
the energy-momentum tensor in Einstein equations should be zero. Therefore, by 
consistence, the above vacuum density must be exactly canceled by a bare cosmological 
constant. If we now derive the vacuum contribution in a curved spacetime, we 
should expect, after the subtraction of the bare cosmological constant, a renormalized, 
curvature-dependent vacuum energy density. For instance, in the case of an expanding, 
spatially isotropic and homogeneous spacetime, filled with vacuum and matter 
components, the renormalized A should be time dependent, being very high at early 
times, but decreasing to zero or nearly zero as the universe expands [2]. 

To have an idea of what a varying cosmological term may be, let us initially 
consider the case of a de Sitter spacetime. It is generally believed that the de Sitter 
cosmological horizon has an associated temperature given by H /2ir, where H = a/a is 
the Hubble parameter [3J1- Therefore, a phenomenological expression for the effective 
vacuum density in this case may be derived by substituting m + H for the energy scale 
m in (|2|), with a subsequent subtraction of Ao « m 4 . In this way we obtain 

A » (m + H) 4 - m 4 . (3) 

In the limit H >> m, this leads to the cutoff-independent result A w H 4 . Since, 
in a de Sitter spacetime, A ss H 2 , this implies that ffisl, that is, this limit necessarily 
describes a de Sitter universe with horizon radius of the order of the Planck length. 
On the other hand, in the limit H << m we have A ss m 3 H. Using again A m H 2 , 
we obtain H m m 3 , or A w m 6 . If, as discussed above, we choose m as the energy 
scale of the QCD vacuum phase transition (of the order of the pion mass), the first 
result is an expression of Dirac's large number coincidence [7] , while the last one gives 
approximately the current observed value of A [8] . 

The above discussion concerns stationary spacetimes, for which H and A are 
truly constants. Nevertheless, it suggests the possibility of a universe evolving from 
an initial, asymptotically de Sitter phase, with A ~ 1, to a final, asymptotically de 
Sitter phase with A ~ m 6 << 1, with A decreasing with time according to ([3]). Before 
verifying such a possibility in the following sections, let us briefly discuss the energy 
conservation in this context. 

There is a common belief about the impossibility of a varying cosmological 
term, because of the Bianchi identities and the covariant conservation of the energy- 
momentum of matter. Indeed, the Bianchi identities G^. M = (G is the Einstein 
tensor) imply, via Einstein's equations, the conservation of the total energy-momentum 
tensor, Tjf.^ = Tjf.„ + g^A iAl = 0, where T is the energy-momentum tensor of matter. 
If one assumes the independent conservation of matter, i.e., = 0, it follows that 
A,u = 0, that is, A is a constant. 

However, if matter is not independently conserved, A may vary with time [9]. In 
the realm of a FRW spacetime, the conservation of the total energy-momentum tensor 
reduces to the form 

p T + 3H(p T +p T ) = 0, (4) 

where px and pt are the total energy density and pressure, respectively. By using 
Pt = Pm + A and pr — Pm — A (p m and p m refer to the corresponding matter 

§ This result was derived by Gibbons and Hawking on the basis of Euclidian methods [3], but the 
positiveness of the de Sitter temperature depends on some appropriate physical interpretation (see, 
for example, \E\ [6]). 
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Figure 1. The Hubble parameter as a function of time (in Planck units) 



quantities), we have 

p m + 3H(p n 



Pm) = "A. 



(5) 



In other words, the vacuum decay is concomitant with a process of matter production, 
in order to preserve the covariant conservation of the total energy. 



2. The early times 

In the spatially flat case, the Friedmann equations give pr — 3H 2 . In the limit of very 
early times we can take H >> m, and ([3]) reduces to A = 3-ff 4 , where the factor 3 
is not important and was chosen for mathematical convenience. Using for matter the 
equation of state of radiation, p rn — p m /3, the conservation equation ([5]) then leads 
to the evolution equation 



H 



2H 2 - 2H 4 = 0. 



For < H < 1, the solution of ^ is given by 

2t = — -tanh -1 il, 
H 



(6) 



(7) 



where t is the cosmological time, and an integration constant was conveniently chosen. 

This solution is plotted in Figure 1, with t and H expressed in Planck units. We 
can see that this universe has no initial singularity, existing since an infinite past, when 
it approaches asymptotically a de Sitter state with H = 1. During an infinitely long 
period we have a quasi-de Sitter universe, with H m A w 1. However, at a given time 
(arbitrarily chosen around t = 0), the expansion undertakes a fast and huge phase 
transition, with H and A decreasing to nearly zero in a few Planck time^jj]. 

From p rn = px — A we can obtain Cl m = 1 — H 2 , where fl m = p m /3H 2 is the 
relative density of matter. Its time variation is plotted in Figure 2. One can see 
that f2 m changes suddenly during the phase transition, from nearly zero, in the quasi- 
de Sitter phase, to nearly 1 at the end of the transition. In this way, the present 
solution has some attributes of an inflationary universe, with an infinitely long period 



| The time at which the transition occurs depends on the integration constant in l(7}. Note, however, 
that it takes place a definite time before the present time. In this sense, we can attribute an age for 
the subsequent universe, which will be determined below. 




Figure 2. The relative energy density of radiation as a function of time 



of inflation ending in a fast transition during which the vacuum decays, releasing a 
huge amount of energy in form of radiation and relativistic matter. After the transition 
we have a radiation-dominated FRW universe, whose subsequent evolution, as we will 
see, is similar to the standard ACDM recipe. 



3. Late times 

In the opposite limit we can take H « m, and ([3]) reduces to A w m 3 H. Let us 
write it as A = oH, with a w m 3 , and let us introduce the equation of state of matter, 
Pm — (7 — l)Pm- From ([5]) we now have 

2H + 3-/H 2 - ajH = 0. (8) 

For H > and p m > 0, one has the solution [TUj 

a(t) = C [exp (cryt/2) - 1]* , (9) 

where C is an integration constant, and a second one was conveniently chosen. 
From it we can derive H (a) and, with the help of p m = 3H 2 — A, we obtain 



Pm — 



A = 



3 (a) 



37/2 



1 



37/2' 



37/2' 



(10) 



(11) 



3.1. The radiation era 

In the radiation phase, doing 7 = 4/3 and taking the limit at << 1, we have 



a w y / 2C 2 at/3, 

4i 2 ' 



Prn 



and 



A = 



3a 4 

o 2 C 2 
3a 2 



a 

2b' 



(12) 
(13) 

(14) 
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The two first results are the same we obtain in the standard model. The third 
one shows that, in this limit, A is sub-dominant compared to p m , and, therefore, the 
matter production can be neglected. This is the reason why radiation is conserved, 
with an energy density scaling with a~ 4 . 

Equations (fT2|) and (fT3j) guarantee that physical processes taking place during the 
radiation phase are not affected by the vacuum decay. For example, the primordial 
nucleosynthesis remains unchanged, since the expansion and reaction rates are the 
same as in the standard context. 

3.2. The matter era 

In the case of a matter fluid dominated by dust, taking 7 = 1 we obtain, from ([9]), 

a(t) = C*[exp(crt/2) - 1] 2/3 . (15) 

In the limit at << 1, it reduces to 

a(f) = C{at/2) 2 ' 3 , (16) 

that is, the scale factor evolves as in a dust-dominated Einstein-de Sitter universe. 
On the other hand, for t — > 00 we have 

a(t) = C*exp(cri/3), (17) 

i.e., our solution tends asymptotically to a de Sitter universe, with H = a/3. 

For the matter and vacuum energy densities, we obtain, from (jlOp and (|11|) . 

a 2 C 3 a 2 C 3 / 2 , 1Q , 

p " l -^r + ^7T' (18) 

A = y + w (19) 

It is not difficult to interpret these results. The first terms in these equations are 
the standard expressions for the scaling of matter and the cosmological term, valid 
if there is no matter production. The first term of (|18p is dominant in the limit 
a/C << 1, when a scales as in (fTl?)) . The first term of (fH?|) . on the other hand, is 
dominant for a — > 00, acting as a genuine cosmological constant. 

The second terms in (|T5)) and (fTO)) are related to the vacuum decay. For large 
times the matter density decreases slower then in the ACDM model, because of the 
matter production. 

3.3. Supernova constraints 

As we have seen, the cosmological scenario presented here is similar, on a qualitative 
level, to the standard scenario of cosmic evolution. Nevertheless, the vacuum decay 
constitutes a substantial difference at late times, and we should verify its consequences 
for the quantitative determination of cosmological parameters like the relative density 
of matter and the universe age, for example. For this purpose, the analysis of the 
redshift-distance relation for supernovas la is of particular importance [TT] . 

From (fT5")) . it is easy to derive the Hubble parameter as a function of the redshift 
z = do/a — 1, where do is the present value of the scale factor. One obtains 

H(z) = H [l - O m + Q m (l + zf' 2 ] , (20) 
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Figure 3. Hubble diagram for 115 supernovae from SNLS Collaboration |11| . 
The curves correspond to Hq — 70 Km.s ^.Mpc , and selected values of f^m- 
For the sake of comparison, the flat ACDM scenario with Q m = 0.27 is also shown. 



where Hq and f2 m are the present values of the Hubble parameter and of the relative 
density of matter, respectively. 

One can see that, as in the ACDM model, we have two parameters to be adjusted 
by fitting the supernova data. In order to do it, we have used (|20|) to fit the data of the 
Supernova Legacy Survey (SNLS) Collaboration |12] , The result is shown in Figure 
3, where we have plotted our theoretical redshift-distance relation for three different 
values of f2 m , together with the theoretical relation predicted by the spatially flat 
ACDM model with fi m = 0.27. In all cases we have used H Q = 70 (Km/s)/Mpc. 

The best fit obtained with ([20]) is given by fl m = 0.32 and H Q = 70 (Km/s)/Mpc, 
with xl = 1-0. With 95% of confidence level, we have 0.27 < O m < 0.37 and 
0.68 < h < 0.72 (h = iJo/lOOKm.s^.Mpc- 1 ). 

With these results we can estimate the universe age in this model. From (|15p we 
can derive an age parameter given by 

_ §ln(n m ) 

Hoto - ~n—T- (21) 

By using the obtained values of Hq and f2 m , one has to ~ 15.7 Gyr, corresponding to 
an age parameter H t — 1.12. 

It is also possible to obtain from (|15p the present deceleration parameter, whose 
best value is q = —0.52. The redshift of transition between the decelerated phase and 
the accelerated one is z = 1.62 [11], showing that we have a decelerated phase long 
enough to permit structure formation. 
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4. Conclusions 

We have described, on the basis of a macroscopic approach to vacuum dynamics, a non- 
singular cosmological scenario in agreement with our general standard view about the 
universe evolution. We have an initially empty, inflationary spacetime, which is driven 
to a radiation-dominated phase through a fast phase transition during which a huge 
amount of energy is released at the expenses of the vacuum decay. The radiation phase 
is indistinguishable from the standard one, and it is followed by a matter-dominated 
era, which evolves to a final de Sitter phase, with vacuum dominating again. 

This scenario is also consistent with a quantitative analysis of the observed 
Hubble diagram for supernovas of high redshifts, leading to cosmological parameters 
in accordance with other - non-cosmological - estimations. In particular, the matter 
density and age parameters are in the intervals imposed, respectively, by dynamical 
estimations of dark matter [13] and globular clusters observations [14] . 

The reader may object that our matter density parameter is above the values 
estimated on the basis of current observations of the cosmic background radiation 
[15] and barion acoustic oscillations [16] - despite the superposition of the intervals 
of allowed values. Nevertheless, it must be emphasized that such estimations are 
dependent on the adopted cosmological model, and should be redone in our case. 
In particular, the production of matter modifies the standard relations between the 
dynamic parameters at the time of last scattering and the present matter density. An 
analysis of such observations in the context of the present model is in progress. 

Another point to be analyzed is the evolution of density perturbations, which 
could be modified by the matter production. A preliminary investigation has shown 
no important difference in the growing of the density contrast in the matter era, while 
in the radiation phase the matter production is irrelevant, as discussed above. 
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